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THE NONLINEAR 2D SUPERCRITICAL INVISCID SHALLOW WATER 

EQUATIONS IN A RECTANGLE 

AIMIN HUANG, MAD ALINA PETCU, AND ROGER TEMAM 


Abstract. In this article we consider the inviscid two-dimensional shallow water equa¬ 
tions in a rectangle. The flow occurs near a stationary solution in the so called supercritical 
regime and we establish short term existence of smooth solutions for the corresponding 
initial and boundary value problem. 


1. Introduction 

Motivated by the study of the inviscid primitive equations, we consider in this article the 
inviscid two-dimensional shallow water equations in a rectangle in the so-called supercritical 
regime. It has been shown that a certain vertical expansion of the inviscid primitive equa¬ 
tions leads to a system of coupled nonlinear equations similar to the inviscid shallow water 
equations; see [RTT08b] and [HT14a]. Hence beside their intrinsic interest, the nonlinear 
shallow water equations can be seen as one mode of the vertical expansion of the primitive 
equations. 

The issue of the boundary conditions to be associated with the primitive or shallow water 
equations has been emphasized as a major problem and limitation for the so-called Local 
Area Models for which weather predictions are sought and simulations are performed within 
a domain for which the boundary has no physical significance, so that there are no phys¬ 
ical laws prescribing the boundary conditions (see [WPT97] and e.g. [RTT08a, RTT08b], 
[CSTT12, SLTT]). The choice of the boundary conditions relies then on mathematical 
considerations (derivation of a well-posed mathematical problem), and on general compu¬ 
tational considerations and physical intuition. The boundary conditions suitable for the 
one-dimensional shallow water equations were derived in an intuitive context in the book 
of Whitham [Whi99] and in [NHF08]; see [PT13] for a rigorous study. For general results 
on boundary value problems for quasilinear hyperbolic system in space dimension one see 
[LY85]; for initial and boundary value problems for hyperbolic equations in smooth domain 
see the thorough book [BS07]. The present article follows the study of the one-dinrensional 
inviscid shallow water equations in [PT13, HPT11] and the study of the linearized shallow 
water equations in [HT14a]. In the study of the linearized inviscid shallow water equations 
in [HT14a] we have shown that five cases can occur depending on the respective values of 
the velocity and the height (not counting the non-generic cases and the symmetries). The 
nonlinear case that we consider in this article relates to what was called the supercritical 
case in [HT14a]; see [HT14b] for the study of a subcritical case. 
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In this article, we consider the inviscid fully nonlinear 2D shallow water equations (SWE) 


( 1 . 1 ) 


u t + uu x + vuy + g(j) x - fv = 0, 

< v t + uv x + vv y + g4> y + fu = 0, 

+ u4> x + V(j)y + 4>(u x + Vy) = 0; 


here JJ = (it, i;, 0)*, (x,y) £ Q := (0, L\) x (0, L 2 ), t £ (0, T), it and v are the two horizontal 
components of the velocity, <p is the height of the water, and g is the gravitational acceler¬ 
ation, / is the Coriolis parameter. The first and second equations (1.1) are derived from 
the equations of conservation of horizontal momentum, and the third one expresses the 
conservation of mass. We consider equations (1.1) for certain values of u,v,4> as described 
below, corresponding to a ’’supercritical” flow and we associate with (1.1), initial conditions 
for u, v, (f> and boundary conditions at x = 0 and y = 0, it, v, (f) vanishing on that part of the 
boundary. 

This article is organized as follows. After this introductory section, we derive in Section 
2 suitable density theorems, density of certain smooth functions in certain function spaces 
of Sobolev type. Section 3 is devoted to the modified (symmetrized) SWE operator for the 
time-independent and the time-dependent cases. It prepares Section 4 in which we deal with 
the linear SWE, linearized around a non-constant time-dependent flow unlike in [HT14a] 
where the background flow is time-independent. In Section 4 we prove the well-posedness of 
the linearized SWE at the price of a loss of derivatives (see Theorem 4.1), and then the well- 
posedness of the linearized SWE in a short time without a loss of derivatives (see Theorem 
4.2). Section 5 considers the fully nonlinear SWE, for which the local well-posedness result 
is obtained. In the Appendices A and B, we collect some useful theorems about semigroup 
and evolution systems, and several classical estimates about functions in Sobolev spaces. 


2. The density theorems 


In this section, we establish general density theorems for certain Sobolev spaces; the 
results supplement and complement those of Section 3 in [HT14a] which we recall when 
needed. These theorems have independent interest, and also will be needed for proving 
later on that —A generates a quasi-contraction semigroup on certain Sobolev spaces, where 
A is the 2D modified SWE operator associated with suitable boundary conditions. 

Throughout this section, let m be a non-negative integer and let A = X(x, y) satisfy 


(2.1) c 0 < X(x,y) < ci, 

where co,ci are positive constants. Furthermore, we say that A = A (x,y) satisfies the 
positive m-condition (m integer > 0) if A satisfies (2.1) and 


VA £ L°°(fi), for m = 0,1, 
A £ H 3Vm (D), for m > 2, 


where a V b = max(a,6). It is easy to see that if mi > m 2 and A satisfies the positive 
mi-condition, then A also satisfies the positive m. 2 -condition. 

We now set for any function 6 = 9(x,y), TO = A (x,y)9 x + 9 y , where A is assumed to 
satisfy the positive m-condition, and introduce the function space 


A 1 m (D) = {9 £ H m {n) : TO = A 0 X + 9 y £ 
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We observe that A™ (D) is a space of local type, that is 

(2.3) If 0 € € C°°(n), then Oip € A 1 m (D). 

This property follows from TfifO) = r ipT9 + + ipy)9, and A 9 (and hence (A fj x + fpy)G) 

is in H m (ft) because of Lemma B.l i). 

We give an equivalent characterization of the space A{”(f2). In the following and through¬ 
out this article, we let d a = d^dy 2 with a = (ai, 0 : 2 ) and set |a| = a\ + Q 2 . We also denote 
by [<9“, /] the commutator [<9“, f]g = d a {fg) — fd a g. 

Proposition 2.1. We assume that A satisfies the positive m-condition. Then 

(2.4) A[ n (Q) = {6 € LI m (D) : € L 2 (D), V |a| = m}. 

Proof. It is clear that (2.4) holds when m = 0. For m = 1, we observe that 

(2.5) Td a 9 = d a (T0) - (d a X)0 x , 

where d a = d x or d y . Then if 6 € H 1 (fl), (d a X)9 X belongs to L 2 (fl) since VA is bounded, 
and by (2.5), Td a 6 belongs to L 2 (fl) if and only if d a {T6) belongs to L 2 (fl)-, (2.4) follows 
for m = 1. 

For m > 2, we observe that 

(2.6) Td a 9 = d a (T9) - [d a , \}0 X , 

holds for all |ct| = m. We note that for 9 € [3“, A ]9 X belongs to L 2 (D) from Lemma 

B.l in) with d = 2, a = A and k = 3 V m. Hence by (2.6), for |a| = m, Td a 9 belongs to 
L 2 (D) if and only if d a (T9) belongs to L 2 (Q), and (2.18) follows as well for m> 2. 

□ 

Now, we need to show that the smooth functions are dense in A^D). Later on we will 
prove more involved density theorems, showing that if u £ A{”(D) vanishes on certain parts 
of dkl, then u can be approximated in A 1 m (D) by smooth functions, vanishing on the same 
parts of the boundary. For the moment, we prove the following: 

Proposition 2.2. C°°(f2) n A 1 m (D) is dense in A 1 m (D). 

Proof. Using a proper covering of Q by sets do, Oi, ■ ■ ■ , On, we consider a partition of unity 
subordinated to this covering, 1 = V ; j ■ Here and again in this section we will use a 

covering of D consisting of Oq, a relatively compact subset of D, and of sets Oi of one of the 
following types: Oi is a ball centered at one of the corners of D, which does not intersect 
the two other sides of D; or Oi is a ball centered on one of the sides of D which does not 
intersect any of the three other sides of D. 

If 9 € A 1 m (D), then 9ipi € A 1 m (D) by (2.3), so that we only need to approximate 9ifi 
by smooth functions. Here the support of t/a is contained in the set Oi, and we start with 
considering the set Oq, relatively compact in ft, then we consider the balls Oi centered on 
the boundary dfl. 

For any function v defined on ft, here and again in the following we denote by v the 
function equal to v in ft and to 0 in M 2 \D. We first consider the case 1 fi = "00 and Oi = Oq 
which is relatively compact in ft. Let p be a mollifier such that p > 0, f p = 1, and p has 
compact support. 
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i) The function v = 9ifo £ A 1 m (H) has compact support in Oq. Since Oq is relatively 
compact in H, then for e small enough, p e *v is supported in H. Using the characterization 
(2.18) for v, the standard modifier theory (see e.g. Appendix C in [Eva98]) shows that for 
e —>• 0: 


(2.7) 


p e *v^-v, in H m (Ll), 

p e * Td a v —> Td a v , in L 2 (D), V |a| = m. 


Since the convolution and the operator T do not commute in the non-constant coefficient 
case, we need the following Friedrichs’ lemma (see [Fri44] or [Hor61, Theorem 3.1]). 

Lemma 2.1. Let U be an open set o/M rf . If Va £ L°°(U) and u € Lf oc (U), then for all 
1 <j<d, 

ad Xj (u * p e ) — ( ad Xj u) * p e —> 0, when e —>• 0, 
in the sense of L 2 convergence on all compact subsets ofU. 


We then continue the proof of Proposition 2.2. Noting that v = 9ifo has compact support 
in 12, we apply Lemma 2.1 with U = Ll,a = A and u = d a v; we obtain 

(2.8) T(Pe * d a v) — p e * Td a v — > 0, V |a| = m, 

in L 2 (D) as e —> 0. Combining (2.7)2 and (2.8), we obtain that as e —> 0, 

(2.9) T(pe * d a v) — > Td a v , in L 2 (D), V |a| = m. 

Therefore, p e *v converges to v in T 1 m (D) by (2.4), (2.7) 1 and (2.9). 

ii) We then consider the case where ^ = if i, and Oi = 0\ which is a ball centered at 
the origin (0,0); the other cases are similar or simpler. Set v = Qif i, and note that v does 
not vanish in general on the boundary dLl of H. In order to extend v to the whole space 
M 2 , we use a well known extension result (see e.g. [Gri85, Theorem 1.4.3.1]). 

Lemma 2.2 (Extension Theorem). Since the boundary dfl of the domain D is Lipschitz 
continuous, there exists a continuous linear operator P = P m from H m (Ll) into H m ( R 2 ) 
such that for all u £ H m (Ll), the restriction of Pu to Ll is u itself, i.e. 

Pu\n = u. 


We denote by v the extension Pv given in Lemma 2.2, and then observe that, for all 
a\ = m, 


( 2 . 10 ) 


d a v = d a v + uf, 
Td a v = Td a v + i/!f , 


where i/f 1 is a measure supported by 0\ \ Q, and uf is a measure supported by 0\ H dLl. 
The two identities in (2.10) together show that 

(2.11) Td a v = Td a v + p a , V |a| = m, 


^In fact, i/f is the sum of a function with support in Oi \ Q and a measure supported by Oi fl dfl, but 
this additional information is not useful to us. 
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where the p a are measures supported by 0\ \ D. Let p be the same modifier as before, but 
now p is compactly supported in {x < 0,y < 0}; then mollifying (2.11) with this p gives 

(2.12) p e * (Td a v) = p e * Td a v + p e * p a , V |ct| = m. 

By the choice of p , p t * p a is supported outside of 0. Hence, restricting (2.12) to 0 implies 
that: 

(2.13) {p e * (Td Q h))| n = Pe * ( Td a v ) ->• Td a v, as e 0, 

in L 2 (D). Applying Lemma 2.1 with U = M 2 , a = A and u = d a v, we obtain that as e —> 0, 

(2.14) T(p e * d a v) — p e * Td a v —> 0, V |a| = m, 
in L 2 (D), which, combined with (2.13), implies that 

(2.15) T{ P e * d a v)\ n -4 Td a v, as e -> 0, 

in L 2 (D). If we set v e = p e * v, then as e —> 0, v e —> v in H m (R 2 ), and 

(2.16) iniL m (0); 

|^T3“(f) e | Q ) —> Td a v , in L 2 (D), V|a| = m, 

which shows that h4 converges to v in A™(D). □ 

Since we have to prove a density theorem involving the boundary values of the functions 
on <9f2, we first need to show that the desired traces at the boundary make sense. We thus 
prove the following trace result. 

Proposition 2.3 (A trace theorem). We assume that A = A (x,y) satisfies the positive 
0-condition. If 0 € 4(D), then the traces of 6 are defined on all of dll, i.e. the traces of 6 
are defined at x = 0,Li, and y = 0,Z/2, and they belong to the respective spaces H y 1 (0,L2) 
and fl'“ 1 (0,Li). Furthermore the trace operators are linear continuous in the corresponding 
spaces, e.g., 9 E 4(^) #|®=o continuous from 4(D) Hy 1 (0,L2). 

Proof. Since 9 E L 2 (D), we see that 9 y belongs to L 2 (0, L 2 )), which implies 

that A 9 X belongs to L 2 (0, Li; H~ l { 0, Lfi)) by observing that T9 = A 9 X + 9 y E L 2 (D). Using 
assumptions (2.1) and (2.2) when m = 0 for A, we obtain 9 X E L 2 (0, H y 1 ( 0, L 2 )), which, 
in combination with 9 E L 2 (0, L\\L 2 (0, Lfi)), shows that 9 E C x ([0, L\\, Hf 1 ^, L 2 )). Hence, 
the traces of 9 are well-defined at x = 0 and L\, and belong to H y ( 0, Lfi). The continuity of 
the corresponding mappings is easy. The proof for the traces at y = 0 and L 2 is similar. □ 

We are now going to introduce density theorems involving the boundary values of the 
functions on OH. Here and throughout this article we denote byFi,r 2 ,T 3 ,r 4 the boundaries 
x = 0, x = Li, y = 0, y = L 2 respectively, and define T to be Ti U T 3 = {x = 0} U {y = 0}. 
We also write 9 |r = 0 as a short notation for 9\ x= q = 9 \ y= 0 = 0, and we introduce the 
function spaces: 

(2.17a) ftF(fi) = {9 E H m (Q) : d a 9\ v = 0, V|a| <m- 1}, 

(2.17b) 4 l (D) = {9 E H m {fil) : T9 = A {x,y)9 x + 9 y E H m {fil), d a 9 | r = 0, V|a| < m}, 
(2.17c) Vr(D) = {9 E C°°(D) : 0 vanishes in a neighborhood of T}. 
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Note that when m = 0, the space is the space L 2 (fil). 

We first have the following characterizations for the space Ap n (n). 

Proposition 2.4. For all integer m > 0, we have 

Tf(H) = {9 £ U?(Q) : TO £ WP(fi), d Q fl| r = 0, V|a| = to}; 

/or all integer m > 1, we /iare 

(2.18) 4^(12) = {0 € WP(fi) : T0 € ^P(fi)}. 

Proof. It is clear that the first statement holds for m = 0, we thus only need to show the 
second statement. By definition of the spaces %p(f2) and Tp ra (12), we see that 

*f*(fi) C {9£ H$(n) : TO £ «P(fi)}. 

In order to prove the converse inclusion, let 9 belongs to the right-hand side of (2.18), then 
it is clear that we only need to show that 9 satisfies the boundary conditions d a 0\ T = 0 for 
all |a| = to. Furthermore, we only need to show that 

(2.19) &?9\ri = 0, and d™9\r 2 = 0, 

since the other boundary conditions involve the derivatives with respect to tangential di¬ 
rections. Since T9 £ T~Cp(Q), on Ti = {x = 0}, we have 

m— 1 

( 2 . 20 ) 0 = d™~\\0 x + 9 y )=Y, d™- k -\\(x,y))d% + d™~ l 9 y + \d™9, 

k =0 

which, together with 9 £ 'Hp(fl), implies that A d™9 = 0. We thus have dfffQ = 0 on Tj. 
Similarly, we can also show that dffQ = 0 on T 2 - We thus completed the proof. 

□ 

As an immediate consequence of Proposition 2.1, we also find the following equivalent 
characterizations of the space Xffiifil). 

Proposition 2.5. 

A'P(n) = {9 £ H m (Q) : Td a 9 £ L 2 (Q), V \a\ = to; d a 9 | r = 0, V \a\ < to} 

= {9 £ : Td a 9 £ L 2 (12), d a 9\ r = 0, V|a| = to}. 

Recall that T = Ti U T 3 = {x = 0} U {y = 0}, and then we state the density theorems: 
Theorem 2.1. Vr(12) fl %p(12) is dense m'Hp(D). 

Theorem 2.2. Suppose that A = A (x,y) satisfies the positive m-condition. Then we have 

Vr(^) H Tp l (H) is dense in 

Remark 2.1. Theorem 2.1 generalizes the classical density results, i.e. thatC°°(kl) is dense 
in and that C“(12) is dense in Hffi(Q), to the functions which vanish on part of the 

boundary <912. 

Theorem 2.2 considers the density of the function space involving functions like 

T9, hence this result is not included in Theorem 2.1. 
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The proof of Theorem 2.1 is similar to or simpler than the proof of Theorem 2.2, we thus 
only prove Theorem 2.2. To prove Theorem 2.2, we proceed similarly as in the proof of 
[HT14a, Theorem 1] 

Proof of Theorem 2.2. Let p(x,y ) be a modifier such that p{x,y) > 0 ,f pdxdy = 1 and p 
has compact support in {0 < \x < y < 2x}. For 9 £ and all a satisfying |a| < m, 

we observe that, 

(2.21) d a 6 = d° r 6 + i/f, TW = TW + i/ 2 “, 

where v f and vf are measures supported by {x = L\} U {y = Li}. Therefore, we have 

(2.22) Td a 9 = Td^9 + p a , V|a|<m, 

where p a = TVf + is a measure also supported by {x = L{\ U {y = L^f. 

We now set 9 e = p e * 9, and mollifying (2.22) with p (see [Hor65]) gives 

(2.23) p e *Td a 9 = p £ *Td°6 + p t * p a , V|a| <m. 

By the choice of p, we have that p e * p a is supported in fl c , and hence restricting (2.23) to 
II implies that 

(2.24) (p e *Td a 9)\ n = (p t *Td«9)\ n ^Td a 9, in L 2 (D) as e -> 0. 

Direct computation shows that 


(2.25) Td a (p e *9) — p e * Td a 9 = x d ^*^ e ) _ Pe * x d ( d * e ) ^ q 


0 , 


dx dx 

where the convergence is in L 2 (D) and achieved by applying Lemma 2.1 with U = R 2 , 
a = A, and u = d a 9. Combining (2.24) and (2.25) yields 

Td a (p e * 0)| n ->• Td a 9, in L 2 (D), as e —>■ 0, 

that is for all |a| < rri. 

(2.26) Td a (9 e |n) -> Td a 9, in L 2 (D), as e 0. 

Similarly, by (2.21), we have for all |a| < m, 

(2.27) d a (9 e \ n ) = (p e *d a 9)\n = (pe*d° r 9)\ n ^d a 9 , in L 2 (D) as e 0, 
where we used the fact that the support of p e * nf is in fl c . In conclusion, there holds 

'9 e \ n ^9, in H m (Q), as e —)• 0, 

Td a (9 e \n) —> Td a 9, in L 2 (Q), as e —> 0, V|a|=m. 


(2.28) 


Finally, 9 t vanishes in a neighborhood of T since the support of 9 t is away from T by 
the choice of p. We thus completed the proof of Theorem 2.2. □ 

Remark 2.2. Looking back carefully at the proof of Theorem 2.2, we see that Theorem 2.2 
is also valid ifT = T 2 UF 4 . Moreover, we say that A (x,y) satisfies the negative in-condition 
if X(x,y) satisfies (2.2) and the following condition: 

(2.1') - ci < \{x,y) < -c 0 , 
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where co,ci are positive constants. Theorem 2.2 is also true if T is Ti U T 4 or T 2 U T 3 , 
and X(x, y) satisfies the negative m-condition provided we choose properly the support of the 
mollifier. 


3. The time dependent shallow water equations operator 

In this section, we aim to study the semigroup property of the (modified) SWE operator 
(see below) with variable coefficients in the supercritical case on the Hilbert space TUff (H) 
(see (2.17a)) with T = Ti UT 3 . We will successively consider the time-independent and the 
time-dependent cases. The linearized SWE operator that we consider reads 

( uu x + vu y + gfi x \ 
uv x + vv y + gfi y J , 

U(f x + V(f y + fi(u x + Vy ) J 

where U = (u,v, ,U = («,?;,(/>)*; we set 

/ u 0 g\ (v 0 0\ 

Si(U) = 0 u 0 , £ 2 (U) = 0 v g ] . 

0 u) \0 4> v) 

Note that £i 2 , £2 admit a symmetrizer So = diag(l, 1, g/(f), i.e. Sq£i, Sq £2 are both symmet¬ 
ric. In order to take advantage of that, we consider the following modified SWE operator: 

(3.2) A°(U)U = £^{U)U X + £$(U)U y , 

where 


£°(U) = S 0 1/2 £i(t/)S 0 - 1/2 = 

^ u 0 


0 u 

0 


1 

0 

u 


(v 0 

° \ 

£$(U) = sl /2 £ 2 (U)Sfi l/2 = 

0 v 



V° \[df> 

V / 


The reason why we choose the form (3.2) will become clear in the next section. 

In the following, we assume that m > 3, the cases when m = 0,1, 2 are similar or simpler. 
Here, we only consider the generic case when U does not vanish, and we first consider the 
time-independent case. We thus assume that U only depends on the space variables x, y 
and that U satisfies the positive (m + l)-condition (see (2.1)-(2.2)) introduced in Section 2, 
i.e. 

(3.3) u, v, <p satisfy the positive (m + l)-condition; 

the reason why we assume one more level of regularity on U will be explained below. As 
we indicated before, we only study the supercritical case, and we thus assume that U also 

2 We sometimes write £1 = £i(U) for the sake of conciseness, etc. 
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satisfies the enhanced supercritical condition: 

(3.4) u 2 -g(f)>cl, v 2 -gcj)>cl, 

where C 2 is a positive constant. 


3.1. Boundary conditions. We aim to determine the boundary conditions which are 
suitable for the system 

(3.5) A°(U)U = £i(U)U x + £ 2 {U)U y = F, 

where F G 'Hp (D) 3 . With assumption (3.4), we see that £® and £ 2 are both positive definite. 
Thus, it is natural to treat either the x- or ^-direction as the time-like direction. Let us 
choose the y-direction, which means that we first need to specify the boundary conditions at 
y = 0 (time-like initial conditions); choosing the x-direction would lead to the same result. 
Multiplying both sides of (3.5) by {£ 2 )~ l gives 

(3.6) Uy + £§(U )- 1 «S 1 °(U)C4 = £° 2 {U)- l F. 


We set kq (U) = \Jg{u 2 + v 2 — gcj))/(j), and we explicitly compute the eigenvalues of (£ 2 ) 1 £ 0 ' 


(3.7) 


uv + 4 >k 0 

Al = ~2 - T’ A2 

v z — gcp 


UV — (j)K o 
v 2 — gcj) 


A3 


u 

v 


Note that all the eigenvalues Ai,A 2 ,A 3 of {£ 3 )^ l £i are positive under assumption (3.4). 
Therefore, from the general hyperbolic theory (see Chapter 4 in [BS07]), it is necessary and 
sufficient to specify the boundary conditions for U at x = 0 in order to solve (3.6) in U. 

In conclusion, in order to solve (3.5) in U, we need to specify the boundary conditions for 
U at x = 0 and y = 0. We then consider the homogeneous case and thus choose to specify 
the boundary conditions for U: 


(3.8) 


U = 0, on T = Ti U T 3 = {x = 0} U {y = 0}. 


As we will see in Lemma 3.2 and Section 5, any sufficiently regular solution U for (3.5) 
and for the nonlinear equations (5.1) will satisfy the following compatibility boundary 
conditions: 


(3.9) 


J d x U = 0, on Ti = {x = 0}, V0 < k < m, 

| dyU = 0, on T 3 = {y = 0}, V0 < k < m, 

which, by differentiating with respect to the tangential direction, is equivalent to 


(3.9') d Q U = 0, on r = T 1 U T 3 , V \a\ < m. 

Hence in the following, we use the compatibility boundary conditions (3.9) rather than the 
boundary conditions (3.8) for the domain of the unbounded operator A defined below. 

We write 77p = 77p(D) 3 for k > 0, in which the functions vanish on T (the part of the 
boundary diV). and we endow the space 77p with the Hilbert scalar product and norm of 


(U,U) n? = J2 (d au ,d a U) L 2 {n) , \\U\\ H? = 

\a\<m 


{(U,U)h?} 1/2 -, 
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we then define the unbounded operator A on W™, by setting AU = A°(U)U, VC/ G P(A) 
and 

v{A) = {u e'H'p : a°(u)u = s\\u)u x + e$(u)u y € n r fi}. 

Note that the compatibility boundary conditions (3.9) are already taken into account 
in the domain T>(A) (see also Propositions 2.4-2.5) since m > 3. We also introduce the 
corresponding smooth function space V(fi) := Vr(H) 3 . Note that V(f2) is dense in Wjf, 
which is a direct consequence of Theorem 2.1. We also have the following results. 

Lemma 3.1. We assume that U satisfies the assumptions (3.3) and (3.4). Then: 

i) V(D) n T>{A) is dense in D(A). 

ii) T>(A) is dense in "Hp . 


Lemma 3.1 is proven below. In order to prove it, we need an equivalent characterization 
of the domain T>(A), which will allow us to use the density results established in Section 2. 
We introduce the notations k, S, V such that 

(v —u 

(3.10) k(U) = Ju 2 + v 2 -g4>, P(U)" 1 = 


V —U —K 


U V 



then direct computations give 


(3.11) 


V t £?V = diag( 
V*£%V = diag( 


UK 


+ \fgj>v UK - \[gcj>i 


u 


2 {it 2 + v 2 )k ’ 2 {u 2 + v 2 )k ’ u 2 + v 2 
vk + \J~g^)u vk — \J~g^>u ft 


=: diag(ai,a 2 ,a 3 ), 
= : diag(&i,& 2 ,b 3 ). 


' 2{u 2 + v 2 )k ’ 2{u 2 + v 2 )k ’ u 2 + v 2 ' 

We then rewrite the modified SWE operator as 
V t A°(U)U = V l £l{V~fi x + V t £^(VE) y 
(3.12) = V t £°VZ x + V t £%VZ y + V t £ ( {V x Z + V t £%V y Z 

= diag(ai, a 2 , a 3 )Z x + diag(&i, b 2 , b 3 )E y + V t £i'P x E + V^^VyZ. 


Direct computations also show that ai,bi,i G {1,2,3} are all positive away from 0, and 
thus both £® and £$ are symmetric and positive definite. Using repeatedly Lemma B.l and 
noting that U belongs to H m+1 (Q), we see that £® , £® , k, V , V^ 1 , at, bi,i G {1, 2,3} belong to 
Furthermore, the last two terms from the right-hand side of (3.12) both belong 
to H m (Q) 3 , and also to Wfit since V = V{U) belongs to H m+1 (Q), and that is the reason 
why we impose one more regularity on U. Therefore, saying that A°(U)U belongs to "Hp 
is equivalent to saying that a^i )X + b^y belongs to 'Hp(Q) for all i G {1,2,3}. Hence, the 
equivalent characterization of the domain T>(A) is 


V(A) = {U = VE 


G TVp, a^i 


T bi^i : j 


G «P(fi),VtG 


{1,2,3}}. 


Proof of Lemma 3.1. We remark that the statement ii) directly follows from i) since V(H) 
is dense in TLff. We thus only need to prove i). 
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Using the new characterization of T>(A) and applying Theorem 2.2 with A = ai/bi for all 
i £ {1,2,3}, we see that each component of 3 can be approximated by smooth functions 
which vanish in a neighborhood ofr = TiUT 3 . Then transforming back to the variable 
U. we obtain that U can also be approximated by smooth functions in V(D). The proof is 
complete. □ 


3.2. Energy estimate for the operator A. In the following, we denote by (•, •) the 
L 2 (D)-scalar product, and observe that if < M with m > 3 (see Remark 3.1 

below for the cases when m = 0,1, 2), then the i7 m -norm of the functions £®{U), £®(U) are 
bounded by some constant C(M). Here and again in this section, the constant C(M ) may 
be different at different places, but it only depends on M. Then for U smooth in T>{A), we 
compute 

(-AU, U) n p = -(d a (£°U x + £$U y ),d a U) 

\a\<m 

= J2 -( e?(&*u) x + £$(&*u) y ,&*u) 

\a\<m 

+ -((l da ’£i} u *’d a U) + ([d a ,£°}U y ,d a U)). 

\a\<m 

Integrating by parts on the first summation at the right-hand side of (3.13) gives 


(3.13) 


(3.14) 




d a U,d a U) 


\a\<m 


x=L i 


fLi 


dy- / {£%d a U,d a U) 


y=L 2 


dx 


x=Q y =0 

+ <(£i* + 3 ly)d a U,d a U) 


Using the compatibility boundary conditions (3.9') and that both £® and £® are positive 
dehnite to dispense with the boundary terms in (3.14), we find that (3.14) is less than 

E|a|<m5(( f i,x + f 2 Jd a t/,<9 a U), which is dominated by + 11^11^)11^11^, 

which is finally bounded by C'(A/)||t/||^ m from the Sobolev embedding H 2 (Q) C 

and m > 3. Applying Lemma B.l iii) with k = m on the commutators from the right-hand 
side of (3.13), we obtain that the second summation at the right-hand side of (3.13) is 
bounded by 


(3.15) 


E c (lk?lb»ll f/ 4ffi.i->ll s “ t/ IL* + irabdNUi-.lla'Uk). 

|a|<m 


which in turn is bounded by C(Af)||[/||^ m . 

Gathering the estimates for (3.14) and (3.15), (3.13) implies that 

(3.16) {—AU, U ) n p < C , 1 (M)||[/||^ f , 
which is 

(3.17) ( AU,U ) h p > -Ci(M)||u||^ m . 
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Thanks to Lemma 3.1 i), we conclude that (3.17) holds for all U in T>(A). 

Remark 3.1. In the cases when m = 0,1,2, we can easily check that the energy estimate 
(3.17) for the operator A is also valid. Indeed, the estimate for the boundary terms is 
the same, and the estimate for the commutators are simpler by direct calculation with the 
assumption that U satisfies the positive m-condition (m = 0,1,2 see (2.1)-(2.2)). 

3.3. The surjectivity of ui + A. We set wo = C\(M), where C\ (M) is the constant 
appearing in (3.17), and we prove the following lemma. 


Lemma 3.2. Let oj be a real number which is greater than ojo. Then if F belongs to Tiff 1 
with m > 3, the equation 

(3.18) £i(U)U x + £$(U)U y + wU = F, 

associated with the following boundary conditions 


(3.19) U = 0, on T = Ti U r 3 = {y = 0} U {x = 0}, 

admits a unique solution U in T>(A). 


Proof. Since £f is nonsingular, multiplying by (£f) 1 on both sides of (3.18) gives 

(3.20) Uy + £$(U)~ 1 £?(U)U X + £%(U)~ l uU = £f{U)~ l F. 

Let us treat again the y-direction as the time-like direction; then (3.20) becomes a 1- 
dimensional hyperbolic system. We observe that the boundary (x-direction only) is a regular 
open subset in R, and that the boundary conditions satisfy the uniform Lopatinskii condition 
(see [BS07, Chapter 9] or [CP82, Chapter 7]). Hence, the general results in [BS07, Chapter 
9] (see also [CP82, Chapter 7]) guarantee the existence and uniqueness of a solution U 
for (3.20) and (3.19). Using (3.20) and the boundary condition (3.19) at y = 0, we can 
conclude by induction that dyU \ v= q = 0 for all 0 < k < m. Similar results also holds for the 
x-direction. Therefore, the solution U also satisfies the compatibility boundary conditions: 

id^U = 0, on Ti = {x = 0}, V0 < k < m, 

\dyU = 0, on r 3 = {y = 0}, V0 <k< m, 


since F belongs to 77p (i.e. it satisfies (3.21) with m replaced by m — 1). It remains to 
show that the solution U actually belongs to Tiff if F belongs to Tiff. 

For all 0 < |a| < m, we deduce from (3.18) that d a U satisfies the following equations 

(3.22) £i(d a U) x + £$(d a U) y + ud a U = d a F - [d a ,£^]U x - [d a ,£$}U y . 

Taking the L 2 (Q) scalar product of each side of (3.22) with d a U and integrating by parts, 
we arrive at 


u 


\d a U\ 


L 2 m + 2 


rL 2 


(£fd a U,d a U) 


(3.23) 


dy+ f ' (£fd a U,d a U) 

x=0 J 0 


y=L 2 


dx 


y =o 


= (, d a F , d a u) + -((£\\ x + £ly)d a U , d a U) 


- ([ d a ,£"]u x ,d a u) - {[d a ,£"]Uy,d a u). 
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The compatibility boundary conditions (3.21) and the fact that £® and £® are both positive 
definite imply that the boundary terms in the left-hand side of (3.23) are nonnegative, 
and thus the left-hand side of (3.23) is larger than w||d“£/||^ 2 (^)- For the right-hand side 
of (3.23), we use the Cauchy-Schwarz inequality to estimate the first term and the same 
arguments as for (3.14)-(3.17) to estimate the last three terms; then summing (3.23) for all 
|a| < m yields: 

+ Ci(M )||^||^m 
11 M 'ht rn + w 0 ||t/|| 

']~L rn ’ 

with ojq being a constant depending only on M. This implies that U belongs to TL p by the 
assumption w > uq. Finally, AU also belongs to "Hp since d a (AU)\r = 0, V|a| <m — 1 and 

AU = £®(U)U X + El {U)U y = F — (jjU. 

Therefore U belongs to T>(A ), and the proof is complete. □ 

3.4. Semigroup. We now set B = uiq + A, with T>(B) = P(A); then (B, T>(B)) is a positive 
operator on ?7p by virtue of (3.17), and u + B (=ui + ujq + A) is surjective for all uj > 0 
thanks to Lemma 3.2. Hence, Theorem A.l (the Hille-Yosida theorem) implies that the 
operator —B generates a contraction semigroup on 77p, and we then obtain the following 
result as a consequence of Theorem A.3 (Bounded Perturbation Theorem I). 

Theorem 3.1. The operator (—A, T>(A)) generates a quasi-contraction semigroup (R(t))t> o 
on TL'y = 7ip(D) 3 satisfying ||-R(t)|| < e w ° f , Vi > 0. 

Remark 3.2. The constant ujq in Theorem 3.1 only depends on the H m -norm of U. 

3.5. Time-dependent modified SWE operator. In this subsection, we consider the 
case where U also depends on the time variable f, and we impose the following assumptions 
on U: 

(3.25a) U (f) satisfies the positive (m + 2)-condition for all t € [0,T], i.e. U(t) belongs to 
H m+2 (Q) [m > 3) and it satisfies the condition (2.1) with co,ci independent of 
t € [0, T]\ 

(3.25b) U belongs to C([0, T]; H m+1 (Q)); 

(3.25c) U satisfies the supercritical condition (3.4) with C 2 independent of t E [0,T]. 

Under these new assumptions, we see that the unbounded operator —A defined in Subsection 
3.1 generates a strongly continuous semigroup with the same arguments as above, once we 
treat the time variable t as a parameter. To be more precise, we define a family of unbounded 
operators A(t)U on the Hilbert space H with A(t)U = A°(U(t))U, VU G D(A(t)\n) and 

V{A{t)\ H ) = {U € H : A°(U(t))U = £^U(t))U x + £l(U(t))U y € H}, 

where H = 77p(= 77p(D) 3 ) and k can be either m -1, mormf 1. 

Since the positive (m + 2)-condition implies the positive (m + 1)- and m-conditions, we 
thus obtain the following result as an immediate consequence of Theorem 3.1. 


(3.24) 


u 


\u\ 


2 

\nv 


< 


< 
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Corollary 3.1. The operators {— A(t)}t generate quasi-contraction semigroups (Rt,i(s)) s >o 


on TL 


m —1 


, (Ms)) o on 74p, and )R>t,‘ 2 (^))s * > o on TL-^ , and they satzsfy 

C(Ml)s Mt(s)\\ < e C(M)s , H^WII < e c ^ s , 


uM«) < e 


for all s > 0, where M\ is the norm of U in C([0,T]; H m 1 (£l)), M is the norm of U in 
C([0, T]; H m (£l)), and M 2 is the norm of U in C([0, T]; H m+1 (Q)). 

Furthermore, with Remark A.l, it is clear that the family {—A(£)} tg [ 0 ,T] Kato-stable 
(see Definition A.5) in all these three spaces 74p~ 1 ,74p an d %p +1 . 


By the definition of these spaces, the two embeddings 

u™ +1 ^ ^ 

are dense and continuous. Using Theorem A.2, we obtain that 74p (resp. ?4p +1 ) is —A(t) 
admissible for all t G [0, T] with respect to "Hp -1 (resp. Wff). That 74 p C V(A(t) 
and 74p +1 C T>(A(t){wp) holds for all t G [0,T] is clear from the definition, and finally, that 
the mapping t i-a —Aft) is continous in the £(?4p, 74p _1 )-norm or £(74p +1 ,74p)-norm 
follows from the second assumption on U (see (3.25b) above). In conclusion, we find the 
following result. 

Lemma 3.3. Assume that U satisfies the assumptions (3.25a)-(3.25c) and m > 3. Then 
the family {—4l(t)}te[o,r] satisfies the Kato-condition (see Definition A.5) with X = 74p _ , 
Y = 74p or X = 74p , Y = R™ +1 . 


4. The linear shallow water system 

In this section, we aim to study the well-posedness of the linear shallow water system 
in certain Sobolev spaces using the evolution semigroups technique. Keeping the notations 
introduced in Section 3, the linear shallow water system reads in compact form 

(4.1) U t + £ 1 (U)U x + £ 2 (U)U y +£(U) = F, 

where t(U) = (—fv,fu, 0)*, and / is the Coriolis parameter. Note that F which does 
not appear in the linearized shallow water system (1.1) is added here for mathematical 
generality and also for the study of the non-homogeneous boundary conditions or for the 
nonlinear case. Observe that the system (4.1) is Friedrichs symmetrizable (see Chapter 1 

in [BS07]) with symmetrizer So = diag(l, 1 ,g/4>), and in order to take advantage of that, 

~ 1 

we make as before a change of variables by setting U = SqU and substitute into (4.1); we 
obtain a new system for U which reads 

(4.2) U t + A°(U)U + B(U)U = SjjF, 
where 

A°(U)U = £i(U)U x + £$(U)U y , 

B(U)U = sl {(Sp) t + £i(U)(SP) x + £ 2 (U)(S~ h -) y )U + £(U). 


(4.3) 
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If we now assume that U satisfies the conditions introduced in Subsection 3.5, then 
the family of operators {—A(f)} tg j 0i T] satisfies the Kato-condition (see Lemma 3.3). If 
we further assume that U belongs to C([0, T\; H m+2 (Q)) D C 1 ([0, T]\H m+1 (D)), then the 
operators {-B(Q}te[o,T] defined by B(t)U = B(U)U are bounded operators on all the three 
spaces B^ 1 iBjf ,B™ +1 by using the estimates in Lemma B.l. Therefore, with Theorem 
A.4 (Bounded Perturbation Theorem II), the family of operators {— Aft) — .B(t)}te[o,T] is a 
Kato-stable family, and furthermore, we have 

Lemma 4.1. The family {— A(t)— -B(Q}te[o,T] satisfies the Kato-condition with X = BTfi~ , 
Y = TLy or X = W£, Y = H™ +1 . 


Combining Theorem A.5 and Lemma 4.1, we obtain an evolution family on By~ 1 and 
another evolution family on B’ff . From the uniqueness in Theorem A.5, we see that these two 
evolution families coincide on B^. Then this unique evolution family satisfies ( E\ ) — (Efi) 
with X = B'fi ^ 1 and satisfies (Efi) — (E5) with Y = TUfi: (see Theorem A.5 and A.6). Using 
Theorem A.6, we obtain that the following system 


^ = -{A{t) + B{t))U{t) + S*F{t), 
U( 0) = Uo 


admits a unique solution U = U(t) £ C([0,T];'Hp ) if Uo £ By andF = F{t) £ C([0, T]; Wfi )• 
Transforming back to the original variables, we obtain the following: 

Theorem 4.1. Let there he given Uq £ W-fi = Hy(LI) 3 and F = F(t) £ C([0, T]; Bp)■ We 
also assume that the U(t ) are given for all t £ [0, T] such that 

(1) U{t) satisfies the positive m+2-condition for all t £ [0,T] , i.e. U(t ) belongs to H m+2 (Ll) 
(m > 3) and it satisfies the condition (2.1) with co,ci independent of t £ [0, T], 

(2) U belongs to C([0, T]; H m+ 2 (n)) n C 1 ([0,Tj; H m+1 (n)), 

(3) U satisfies the supercritical condition (3.4) with C 2 independent of t £ [0,T], 

Then the system (4.1) associated with the initial condition U{ 0) = Uq has a unique solution 
U = U{t) which belongs to C([0,T];ftJ?) n C 1 ([0,T]; ftJT 1 ). 

Remark 4.1. Using the system (4.1), the solution U in Theorem fil satisfies the compat¬ 
ibility boundary conditions (3.9) by the same argument as in Lemma 3.2. 

Remark 4.2. Notice that ifUo £ 3 , and U,F £ C°°([0, T]; then the solution 

U provided by Theorem 4.1 belongs to C([0,T]; H k (Ll) 3 ) for all k> 0, which implies that U 
belongs to C([0, T]\ C 00 (D)), and then by using the system (4.1), we conclude by induction 
that U actually belongs to C°°([0, T];C°°(f2)). 

We lost two space derivatives from U to the solution U for the linear system (4.1) in 
Theorem 4.1, which is not sufficient for us to study the nonlinear case. In order to gain 
these two derivatives back, we need some additional a priori estimates. 


3 For the domain SI = (0, Li) x (0, L 2 ), we have that C°°(S1) = (~lj *2 0 H k (£i), see [Gri85, Chapter 1], 
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4.1. A priori estimates. With Remark 4.2, we assume that D, Do, F,U are smooth func¬ 
tions satisfying the following system 


(4.5) 


, u t + s 1 (u)u x + e 2 (u)Uy + e(u) = F, 
< u(0 ) = Do, 

U |r = 0, 


with Do, D( t),F(t ) satisfying the compatibility boundary conditions (3.9) (i.e. (3.9')) for all 
t € [0, T], and r = {x = 0} U {y = 0}. In addition, D is positive away from 0 and satisfies 
the supercritical condition (3.4), i.e. 


(4.6) 


co < D < ci, 

u 2 — gcj) > c|, v 2 — gcj) > c. 


• 2 > 


where co,ci,C 2 are positive constants. We will first derive L 2 a priori estimates for the 
linear system (4.5) and then extend the D 2 -estimates to H m -estimates with m > 3. For 
the sake of simplicity, we write fix = H x [0, T], and = D°°(0, T; H k (£l)) for all 

k = 1, • • • , m and L°°(L 2 ) = L°°(0, T ; D 2 (fl)). We assume that 

(4.7) ll^llz, 00 ^™) — ^ M. 

Multiplying (4.5)i by So and taking the scalar product in D 2 (H) with D gives 

(4.8) (S 0 U t , U) + (S 0 £\U X + S 0 £ 2 U y , D) + (S 0 £(U), U) = {S 0 F , D). 

We now calculate 


(4.9) 


(S 0 Dt, D) = ^<S 0 D,D) - ^(S 0 yD,D), 


{So£(U), U) = 0; 

and, using integration by parts, we find that 

r L 2 


{Sq£\U x + S 0 £ 2 U y , D) = l r (S 0 S,D, D) 'dy + i H (S 0 £ 2 U , D) 

^ 40 a:=0 2 •'0 


1 /' L1 


y=L 2 


dx 


y =o 


(4.10) 


-i(((SoS 1 ), + (SoS 2 ) J/ )D,D) 
>-i(((S 0 Si) x + (S 0 S 2 ) y )D,D), 


where the last inequality results from the boundary conditions (4.5)3 and the fact that 
So Si, Sq £ 2 are both positive definite. 

Finally, we obtain the following inequality by gathering the calculations (4.9)-(4.10) : 
(4.11) ^<S 0 D,D> ^ ( S o,tU,U) + (((S 0 £i) x + (S 0 £ 2 ) y )U,U) + 2 (S 0 F,U). 
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11 ~ 112 

We set Io(t) = (SqU, U) = an d then the first two terms in the right-hand side 

of (4.11) are bounded by 

\\{S 0 ,t + (S 0 £i) x + (Soh^S^W (S 0 U,U) 

(4,12) ^ ^ ^ 

— C(M Loo{nTy ||t4|| L oo(Q T )) ||^/|Loo(Q T )l II ^11 L 00 (Q t )^ 0 ^)’ 

which is dominated by C'(||[/t|| i<X) . ii - 2 ^, ||t/|| Loo ^ 3 p/o(i) by using the Sobolev embedding 

H 2 (£l) C L°°(Q). Using the Cauchy-Schwarz inequality, we estimate the last term in the 
right-hand side of (4.11): 

{4 2 ( s „ f , u ) < n4fiii a(n) +iis 0 U||y (n) 

Combining with (4.12) and (4.13), (4.11) implies that 

(4 14) ~dt~ + ^ , (||^'IL oo (J? 3 ))|l^ 1 WlL 2 (fi) 

< r i(U)(l 0 (t) + \\F(t)\\ 2 L2m ), 

where the constant r\{U) = C , (||^t|| i0O (jj 2 p ll^ r IL°o(H 3 )) + ^ on ly depends increasingly on 
ll&ILoo^). ||E/|| i<x>(if3) . We observe that 

ri(U) = ri{\\U t \\ LOO(Hm _ ly \\U\\ Leo(Hm) )<r 1 (M,M) = r 1 (M), 

with m > 3 by the assumption (4.7), and we write r± = r\(M) > 1 for the sake of simplicity. 
Using Gronwall’s lemma for (4.14), we obtain 


(4.15) 


Noticing that 


/o(i) < e rit (/ 0 (0)+n / ||E( S )||^ 2(n) d S ) 


< 


e rtt ri-(\\U4l Hn) +j‘\\F(s)\\l %a) ds). 


Io(t) = + y2 + ^<J> 2 )dxdy > min(l, 5 /ci)||C/(t)|| 42(n} , 

and setting r 2 = l/min(l, g/ci), (4.15) implies that 
(4.16) ||E/(t)fe (n) < e’'“-'ir 2 (||t/„|| 2 L2(n) + J‘ ||F( S )|| 2 1(n) *). 

Taking the L°°-norm of (4.16) over [0, T] immediately gives 

(4.U) HL,„, < C„(M,T)(||(7„|| 2 12(n) +T||F|| 2 „ (13) ), 

where Cq(M. T) = e ri (^) T r 1 (M)r 2 , only depends on the bound of the L°°(H m )-norm of U 
and the L°° norm of Ut. 


18 


A. HUANG, M. PETCU, AND R. TEMAM 


We now turn to extending the L 2 -estimate (4.17) to //"‘-estimate. Applying d a = d^dy 2 
with |a| = ai + a 2 < m to (4.5) and recalling that U satisfies the compatibility boundary 
conditions (3.9 / ), we obtain that d a U satisfies the following equations 


(4.18) 


' (d a U) t + £i(t/)(0°l7)* + S 2 (U){d a U) y + l(d a U) = F a , 
< d a U{ 0) = d a U 0 , 
d a U | r = 0, 


where F a = d a F — [d a , £\}U X — [d a , £ 2 ]U y . Observing that (4.18) has the same form as (4.5), 
therefore proceeding exactly as for (4.17), we find 

(4.19) !|0”t/|l!„ (ia) < C„(M,r)(||8“(7„|| 2 [a(n) + r||F„||l„ (ia) ), 

where Co(M, T) is the same as in (4.17). 

We now need to estimate F a . Lemma B.l iii) with k = m on the commutators in F a 
gives 

+c(||8i(&(*))||U (n) ||cyt)||J, w _ I(n) 

(4 - 20) + ll £ 2(^( t ))||ff~(fi)ll f/ »( t )|llfW- 1( n)) 

< ||a“^(t)||L (n) + c'(xf)||c/( i )||^ l(fi) , 

where C{M) only depends on M-the bound of the norm of U. 

Summing (4.19) for all |a| < m and using the estimates (4.20) for F a , we finally arrive 
at 

< 4 - 21 ) < C„(M,T)(||f/ 0 || 2 Hm(n) +r||F|| 2 1 „ (H „ ) +TC(M)||C/|| 2 i „ (H „ ) ), 

where the constants C(M ) may be different at different places, but they enjoy the same 
property, i.e. they only depend on the bound of the L°°(H m )- norm of U in an increasing 
way. 

We choose T small enough so that Co(M,T)TC(M) < 1/2; with this choice of T, we are 
able to absorb the term ||t/|| Loo ^ m ^ in the right-hand side of (4.21) and we find that 

(4.22) ||C||L (Hm) < 2C„(M,T)(||f/„|| 2 ifm(n) +T||F|| 2 „ (h „ ) ), 

where Cq(M,T ) is the same as in (4.17). We emphasize the fact that the choice of T only 
depends on the bound M of the L°°(H m )- norm of U and the L°° norm of Ut- 
Finally, we estimate the L°°(i7 m-1 )-norm of Ut- We write (4.5)i as 

(4.23) U t = F- £(U) - £i(U)U x - £ 2 (U)U y: 

We first take H m ~ ] (D)-norm of (4.23) and use Lemma B.l i) with s = m — 1 and d = 2 
to estimate the last two terms in the right-hand side of (4.23); then we take L°°-norm over 
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[0, T], and we find 

(4.24) + y||^|| J^OO — + ( ^||^2(U)|| Loo ^ m _ 1 ^||Uy|| ioo( . i?m _ 1 j 

— II^IL^Off™- 1 ) 4“ ^(|l ’ ^)|I^IL°°(.ff™)> 

where / is the Coriolis parameter. The inequality (4.24) shows that 

(4.25) II^IIl^RT™- 1 ) — ^II'^IIl 00 ^’"- 1 ) 4" 2C'i(M)||l/|| i00 (^ i - m ), 

where C\(M) only depends on M- the bound of U in L°°(H m ). We also obtain the following 
L°°(L 2 )-estimate 


(4.26) INlL ( y) £ 2||F|| 2 t „ (i2) + 2C 1 (Af)||(/||h ( „. r 

4.2. Improved regularity. With the LI m -estimates (4.22) and (4.25) at hand, we are now 
able to gain back the derivatives lost in Theorem 4.1 by shrinking down the time T, and 
we prove the following theorem. 


Theorem 4.2. Let there he given 

U 0 £H m (n), F,U € L°°(0,T;H m (n)), U t e L°°(0, T; ^(fl)), 

and furthermore we also assume that Uo,F(t) satisfy the compatibility boundary conditions 
(3.9) for t £ [0, T\, and that U satisfies (4.6) and (4.7). Then there exists T > 0 small 
enough depending only on the bound of the L°°(H m )-norm of U and the L°°(H m ~ 1 )-norm 
of Ut such that the system (4.1) associated with the initial condition U{ 0) = Uq and the 
homogeneous boundary conditions (3.8) has a unique solution U = U(t ) such that 

U € L°°(0,T-,H m (n)), U t € L°°(0, T; L7 m_1 (fl)), 

and the solution U satisfies the compatibility and boundary conditions (3.9) and the esti¬ 
mates (4.22) and (4.25). 


Proof. Let p(x,y), cr(t) be modifiers such that p(x,y),a(t ) > 0, f pdxdy = f adt = 1 and 
p has compact support in {0 < \x < y < 2x}. For a function w defined on D, (p e * w)\q 
stands for the restriction to Ll of p t * ui, where w is the extension of w by 0 outside 12, and 
similar notations are also used for the functions defined in Qt, or the vector functions (with 
the notation applied to each component of the vector functions). We then set 


Uq = (p e * Uo)\n, F e = ((pa) £ *F)\ n , U e = ((pa) e *U )| n . 


Standard modifier theory shows that Uq,F € ,U e converge to Uo,F,U respectively as e 
in the corresponding spaces 4 . Hence for e small enough, we can assume that 


Ik, 




< 2 l|t/o|| 


1^ IL°°(H m ) — 


2 ||i? e || 2 < 2 |IeII 2 

H m ( a)’ II II — II II L 2 (H m )’ 

II Ut < 2AI. 


4 See the details in a related situation in [HT14b] 
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In addition, with the choice of p, we have that the support of Uq is away from r = {x = 
0}U{y = 0}, and so is the support of F e (t) for all t € [0, T], Therefore, Uq, F e (t) also satisfy 
(3.9) for all t € [0, T]. Then using Theorem 4.1 and Remarks 4.1-4.2, there exists a smooth 
solution U e for system (4.5) with Uq,F,U replaced by Uq, F e ,U £ , and U e also satisfies the 
compatibility boundary conditions (3.9). For T > 0 small enough only depending on the 
bound of the norm of U and the L°° (i/ m ^ 1 )-norrn of Ut, then the a priori estimates 

(4.22) gives that 

MlUm £ 2C„(2M,D(||^t m(n) + r||F-|| 2 i „ (ffm) ) 

< 4C„(2M,T)(||f/„|| 2 H . (n) +rj|F|| 2 1 „ w .. ) ). 

where 2 M is the bound of the norm of U e and the norm of Uf. 

The inequality (4.27) gives a uniform bound on the sequence {U e }, which implies that 
there exists a subsequence of {U e } converging weak-star in L°°(0, T; The next 

point is to prove that the sequence {U 6 } is Cauchy in L°°(0, T; L 2 (H)). For that purpose, 
we write 


(4.28) W e = U e -U e ', 

and subtracting the corresponding equations of form (4.5) satisfied by U e and U e ', we obtain 


(4.29) 


' W t e + £i{U e )WZ + £ 2 (U e )W^ + £(W e ) = F, 
< IU e (0) = U^Q - u^, 

W e \ r = 0, 


where 


F = F e - F e ' + (5i (U e ’) - £i{U e ))ui + {£ 2 {U e ') ~ £ 2 (U e ))t //. 


Noticing that (4.29) has the same form as (4.5), therefore proceeding exactly as for (4.17), 
we obtain 


(4.30) ||W'<|| 2 „ (iS) <C 0 (2M,r)(||f/J-f/J|| 2 1(n) + r||F|| 2 „ (i2) ). 
Using the explicit expressions for £\ and £ 2 , direct computation shows that 

(4.31) \\£ 2 {W')-£ 2 (U 
Therefore, combining the estimates in (4.30) and (4.31), we obtain 


2 < 311 f7 e — u e '^ 2 

l°°(l 2 ) - \\ u u 


L°°(L 2 )' 


(4.32) 


\U e - U 


:'||2 


f 'l\2 

I L°°(L 2 ) 

f'||2 

lL°°(n T )^ 




< c„(2M,r)(||t/5 -1/5 || 2 1(n) + r!lF< - f< 

+ 3T||t/*-t/*'|| 2 „ ( y ) (||Cji|U (fiT) 

which, with the use of the Sobolev embedding H 2 (£l) C and noting that m > 3, is 

furthermore bounded by 
(4.33) 


C 0 (2M,T)(||C/ 0 e -C/ 0 e, || 2 2(n) +T| 


F e - F ] 


\ \2 




\L°° 


’ 112 

L°°(L 2 )1 


Since {Uq} is a Cauchy sequence in L 2 (£l), and {F e },{U e } are Cauchy sequences in 
L°°(L 2 ), and {U e } is uniformly bounded in the L°°(77 ,n )-norm by (4.27), we obtain from 

(4.32)-(4.33) that {U e } is also a Cauchy sequence in L°°(0, T; L 2 (fl)). Hence by L 2 — H m 



THE 2D SUPERCRITICAL SHALLOW WATER EQUATIONS 


21 


interpolation, the sequence {U e } converges strongly in L°°(0, T; iL m_1 (fl)) to a function U 
which belongs to L°°(0, T; H m (El)). Using Proposition 2.3, we obtain that U satisfies the 
compatibility boundary conditions (3.9) since U e satisfies (3.9). 

The a priori estimates (4.25) give a uniform bound on the sequence {[//}, i.e. 


(4.34) 


m 


|2 

\L° 


< 9 F e 
'( h ™- 1 ) — z lr 


12 
\L° 


, (ffm _ 1) + 2C' 1 (2M)||[P| 

<4||F||^ (ffm _ 1) + 2C 1 (2M)||[/ f || 


\\ 

2 


since the sequence {U e } are uniformly bounded in the L°°(H m )- norm. 

Proceeding exactly as for (4.26), we obtain that 

(4.35) || Ut - Uf\ |L (I1) < 2||F||^ (i2) + 2C 1 (2M)||tf - U‘'f L „ { „ ly 

which implies that {Uf} is also a Cauchy sequence in L°°(L 2 ) by using the above estimates 
for F and noting that {t/ e } is Cauchy in L°°(H m ~ l ) with m > 3. Therefore, by L 2 — R m_1 
interpolation, we obtain that {Uf\ converges strongly in L°°(H m ~ 2 ) to a function V which 
belongs to L°°(P m “ 1 ). 

Now passing to the limit, we obtain that U solves (4.5), and Ut = V at least in the sense 
of distributions. Finally, proceeding exactly as in Subsection 4.1, we see that the solution U 
satisfies the estimates (4.22) and (4.25); the uniqueness directly follows from the estimate 
(4.22). We thus completed the proof. □ 

Remark 4.3 (Non-homogeneous boundary conditions). Using Remark 9.1 in [HT14a], the 
existence of a solution for the linear system (4.1) associated with non-homogeneous boundary 
conditions can be obtained, we omit the details here; see [RTT08b] for a similar situation. 


5. The fully nonlinear shallow water system 

In this section, we aim to investigate the well-posedness for Eqs. (1.1) associated with 
suitable initial conditions and homogeneous boundary conditions, and we will make a remark 
about the case of non-homogeneous boundary conditions. Keeping the notations introduced 
in Section 3 and 4, the fully nonlinear shallow water system reads in compact form 

(5.1) U t + £ 1 {U)U x + £ 2 (U)U y +£(U) = 0. 

5.1. Stationary solution. We want to study system (5.1) near a stationary solution, and 
we start by constructing such a stationary solution (u,v,(j>) = (u s ,v s ,cj> s ). These functions 
are independent of time and satisfy 

(5.2) £ 1 (U)U x +£ 2 {U)U y + e(U) = 0. 

The existence of the general stationary solution U s to (5.2), which satisfies the supercrit¬ 
ical condition, is a 1-dimensional hyperbolic problem if we multiply by Elf 1 on both sides 
of (5.2) and treat the y-direction as the time-like direction as we already did in Lemma 

3.2. The general results in [BS07, Chapter 11] guarantee the existence of the stationary 
solution U s if we specify suitable initial (y-direction) and boundary (x-direction) conditions. 
Actually, Subsection 2.1 in [HPT11] provides an y-independent stationary solution to (5.2) 
satisfying the supercritical condition. But in what follows, we think of U s in a general form 
(i.e. U s depends both x and y). 
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(5.3) 


(5.4) 


We thus choose our stationary solution U s = (u s , v s , ^ s )* of (5.2) satisfying a strong form 
of the supercritical condition, i.e. 

2Co ^ U s ,V s ,(j) s 'i' 2^1) 

u 2 - g(j) s > 2c|, v 2 - gc/) s > 2c|, 

where Co, ci, C 2 are given, positive constants which will play the same role as those in (4.6). 

We set U = U s + U . Note that if we choose 5 sufficiently small so that if \U\ <5 (i.e. 
|w|, |u|, |<ji| < <5), then U satisfies relations similar to (5.3), that is: 

c 0 < u + u s ,v + v s ,4> + 4> s < Cl, 

% (u + u s ) 2 - g(4> + 4> s ) > c|, (v + v s ) 2 - g(cj) + 4> s ) > c|. 

The relations (5.4) will guarantee that we remain in the supercritical case. 

We then substitute U = U s + U into (5.1); we obtain a new system for U, and dropping 
the tildes, our new system reads: 

(5.5) U t + £i (U + U S )U X + £ 2 (U + U s )Uy + £(U) = F u , 

where 

F u = —£i(U + Us)U s , x - £ 2 {U + U s )U s , y - £{U S ) 

UU S , X + VU S: y 
UV S , X + VV S) y 

T V(f) St y T (j)(u S:X + Vs,y)/ 

In order to have the last equality, we use the fact that U s is a stationary solution satisfying 
(5.2). We supplement (5.5) with the following initial and homogeneous boundary conditions: 

(5.7) I.C. 17(0) = Ho, B.C. U\ x=0 , y= 0 = 0. 

Observe that we can rewrite (5.5) as 

(5.8) U x = £i (U + Us)~\F u -U t - £ 2 {U + U s )U y - £(U)); 

at x = 0, we immediately see that U x \x=o = 0. Applying d x with l<fc<m — lto (5.8), 
we can conclude by induction that d x U | x= o = 0, for all k = 0, • • • , m. Similarly, we also 
obtain that d y U\ y= o = 0 for all k = ().■■ ■ , m at y = 0. Therefore, if U satisfies (5.5) and 
(5.7), then U also satisfies the compatibility boundary conditions (3.9), i.e. 

d x U = 0, on Ti = {x = 0}, V0 < k < m, 


(5.6) 


(5.9) 


d y U = 0, on T 3 = {y = 0}, V0 < k < m. 


5.2. Nonlinear shallow water system. In order to be able to solve the nonlinear system 
(5.5)-(5.7), we require the initial and boundary conditions to be compatible. We thus 
assume that Uq satisfies the compatibility boundary conditions (5.9) (i.e. (3.9)). We are 
now on the stage to prove the following result. 

Theorem 5.1. Let there be given the stationary solution U s £ F[ m+1 (Q) with m > 3, and 
two positive constants Mq,M such that 


Mq, M > 0, 


Mo £ (0,-], M — y Mq + 11C/ s 

Dra 


I 
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where v m denotes the norm of the Sobolev embedding H m ( fl) C L°°(Q). We are also given 
the initial condition JJq £ which satisfies (5.9) and 

(5T°) \\Uo\\ H m m < mm^Mo, 2C {U S ) + 4C , 1 (M)C' 0 (M, 1) ’ 4C 0 (M, 1)) ’ 

where Cq(M, 1) = e ri ^ M \\{M)r 2 (resp. C\{M)) is the constant appearing in (4.17) (resp. 
(4.25)j, and C(U S ) only depends on the bound of the H m+l (fiV)-norm of U s (see (5.18) 
below). 

Then there exists T > 0 only depending on the initial data U$ and the stationary solution 
U s such that the system (5.5)-(5.7) admits a unique solution U satisfying 

U £ L°°(0,T-,H m (n)), Ut £ L°°(0, T; 77 m_1 (0)). 


Proof. As a preliminary, we choose T such that T < 1. 

Considering the compatibility boundary conditions (5.9), the resolution of the nonlinear 
system (5.5)-(5.7) will be done using the following iterative scheme 


(5.11) 


'U t k+1 + £i{U k + U s )U k+l + £ 2 {U k + U s )U k+l + l(U k+1 ) = F u \ 
< U k+1 ( 0) = U 0 , 
k U fc+1 | r = 0. 


We initiate our iteration scheme by setting U° = Uq. and then construct the approximate 
solutions U k by induction. 

If we assume that 

(5.12) ||U fc ||^oo ( 0 iT;fl -m (n)) < M 0 , \\U k \\ 2 LOO{QT . H m-i m < M 0 , U k satisfies (5.9), 

then we have that F uk also satisfies (5.9). Using the Sobolev embedding C L°°(D), 

the L°°-norm of U k is controlled by <5, which shows that U k satisfies the supercritical 
condition (5.4), i.e. 


co < u k + u s ,v k + v s , f k + fi s < Cl, 

(u k + u s ) 2 - g(cj) k + fi s ) > c|, ( v k + v s ) 2 - g((f k + fi s ) > c$; 


furthermore, we have 

|| U k + C7 a || Loo(frm) < y^o + || U s \\ Hm(Q) < M, ||(U fc + U s ) t ILoo^-r) < V // Mo < M. 

Therefore, for T > 0 small enough only depending on M. applying Theorem 4.2 to (5.11) 
with U = U k + U S ,F = F l,k , gives a solution U k+1 which satisfies (5.9), and that 

(5.13) 




and 

(5.14) 


K +1 \\ 2 




H 171 (Cl) 1 II \ \ L°°(H m ) 

< (using Lemma B.l i) with s = m and d = 2 for F u ) 

<2C’„(M,i)(||t/„|| 2 ifm(n) + rc||(7,|| 2 if „ + . (n) ||c/ t ||L (H „ ) ), 
t/ ‘" 2 ,i,™-, ) +C 1 (M)||(7‘ +1 ||L ( „™ ) . 


< 2 T 
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where Cq(M,T) (resp. C\(M )) is the constant appearing in (4.17) (resp. (4.25)), and M is 
the bound of the L°°(H m )-norm of U (= U k + U s ) and the L°°(H r 
It follows from the explicit form (5.6) of F u that 


Tm 1 )-norm of Ut (= Ut ) 


(5.18) 


\u k+1 \r 

\ U t II L° 


(5.15) ( F u ) t = F u K 

In particular, ( F lJk ) t = F u *. Using also that U k \t=o = Uq, we obtain with the mean value 
theorem that 

(5.16) |i^(f)| < | F u \t) -F Uo | + \F U °\ < t\F u t(t')\ + \F Uo \, 
for all t € [0, T] and for some t' € (0, t), which implies that 

\\F uk \\ 2 < 2T 2 llE t/tfc || 2 +‘^lli ?c/ °ll 2 

(5.17) < (using Lemma B.l i ) with s = m — 1 and d = 2) 

< 2r 2 c||c/»||U (!! )IIUIIi„ ( „„-, ) + 2C||c/,||U (n )U(i|lU- 1 ( n )' 

Gathering the estimates (5.13)-(5.17), and using the assumption (5.12), we finally arrive at 

\U k+1 \\l x{Hm) < 2C 0 (M,l)\\U 0 \\ 2 Hmm + TC 0 (M,l)M 0 C(U s ), 

< {C(U a ) + 2C 1 (M)C 0 (M, 1))||C7o|| 

+TMqC(U s ) (T + 2Ci(M)C 0 (M, 1)), 

where C{XJ S ) only depends on the bound of the i7 m+1 (Q)-norm of U s . Note that the first 
two terms in the right-hand side of (5.18) are less than Mo/2 by the assumption (5.10), and 
both the second terms in the right hand side of (5.18) approach 0 when T —>• 0; we thus 
can choose T small enough again such that 

(5.19) \\U k+1 \\ 2 Lao(Hm) < M 0 , K +1 || 2 00(Hm _ 1) < M 0 , U k+1 satisfies (5.9). 

Now U k+l also satisfies (5.12); hence we can continue our construction. Let us emphasize 
that the choice of T only depends on Mo, M, U s and is independent of k, therefore our 
iteration scheme can be conducted for all k. and we can construct the sequence {U k } as 
long as the starting point U° satisfies (5.12), which holds true by the assumption (5.10). 

We now have an uniformly bounded sequence {U k } at hand, and the next point is to 
show that the sequence {U k } is Cauchy in L°°(0, T ; L 2 (H)), which is almost achieved in the 
proof of Theorem 4.2. Let us write W k+1 = U k+1 — U k ; with W e ,U e ,U e ',U e ,U e ', F e , F e ' 
replaced by W k+1 , U k+1 , U k , U k + U s , U k ~ l + U s . F uk , F uk 1 in (4.29), proceeding exactly 
as for (4.32)-(4.33) and noticing that IU fc+1 (0) = 0, we obtain 


| U k+1 - U k f Lao{L2) < C 0 (M,T)T(\\F uk - 


(5.20) 


?U h ~ 1 \ 12 


I L°°{L 2 ) 


+ 6 \\U 


k ||2 


\u k -u k 




\L°°{H m ) IR ^ llL°°(L 2 )l 

Using the explicit expression (5.6) for F u and the Sobolev embedding H m {VL) C L°°(H), 
we estimate 
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Combining (5.20) and (5.21) and using the uniform boundedness of the L°°(0, T; iL m (D))- 
norrn of U k , we obtain 

(5.22) ||(7‘+‘ - E/‘||L (1J) < C 0 (M,T)T(C(\\U,\\ HW , {n) ) + 6M„)||f/‘ - 

Upon reducing T again, we can assume that 
(5-23) C 0 (M,T)T(C(\\U s \\ Hrn{n) ) + 6M 0 ) < i; 

then the inequality (5.22) implies that 

(5.24) ||t/‘+‘ - U k \\ L „ m < I||£7* - <■■■< (l)iC 1 - C»L„ (i2) . 

Therefore {U k } is Cauchy in L°°(L 2 ); let U be the limit of this sequence. Note also that 
{U k } is uniformly bounded in L°°(H m ), so by L 2 — H m interpolation, the sequence {U k } 
converges strongly in L °°to U G Similarly as for (4.35), we can obtain 

that {U k } converges strongly in L°°(H m ~ 2 ) to a function V € L°°(H rn ~ 1 ). 

Now passing to the limit in (5.11), we obtain that U solves (5.5) and Ut = V, and that U 
belongs to L°°(0, T; H m (Q)) and Ut belongs to L°°(0, T ; H m ~ l {p,)). The uniqueness directly 
follows from (5.22). This completes the proof. □ 


Remark 5.1 (Non-homogeneous boundary conditions). With Remark 4-3, the existence 
of a solution for the iterative scheme (5.11) associated with non-homogeneous boundary 
conditions can be obtained, and by passing to the limit, the nonlinear system (5.5) associated 
with non-homogeneous boundary conditions admits a unique solution; we omit the details 
here. 


Remark 5.2. After completing this article, we found that we can also use a finite difference 
method to prove the existence and uniqueness of the fully nonlinear SWE (i.e. Theorem 5.1) 
by observing that we have the energy estimates (3.24) for the corresponding boundary value 
problem (although slightly different), which is the one we only need for the finite difference 
method. We omit the details here. 

The finite difference method has the advantages that we do not need the density theorems 
in Section 2 and the evolution semigroup technique. However, the evolution semigroup 
technique has its own advantage that it tells us how we lost two space derivatives for the 
well-posedness of the linear SWE (see Theorem 4-l)> an d that explains why we only have 
local well-posedness for the fully nonlinear SWE in some sense. 


Appendix A. Preliminary results about semigroups and evolution families 

This appendix collects some basic facts on the semigroups and evolution families and 
the characterization of their generators. The main references are the classical books by K. 
Yosida [Yos80], and by E. Hille and R.S. Phillips [HP74], and by A. Pazy [Paz83] and the 
book by K.-J. Engel and R. Nagel [EN00]. 

Definition A.l. A family (S(t))t >o of bounded linear operators on a Banach space X is 
called a strongly continuous (one-parameter) semigroup (or Co-semigroup) if it satisfies 

i) 5(0) = I , S(t + s) = S(t)S(s ) for all t,s > 0; 

ii) : t i-a £ x (t) := S(t)x is continuous from M + into X for every x € X. 
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Proposition A.l. For every strongly continuous semigroup (S(t)) t > o, there exist constants 
tv € M and M > 1 such that 

(A.l) j|S(t)|| < Me wt 

for all t > 0. 

Definition A.2. A strongly continuous semigroup is called quasi-contraction if we can take 
M = 1 in (A.l), and called bounded if tv = 0, and called contraction if tv = 0 and M = 1 
is possible. 

Definition A.3. The generator A : T>(A) C X > X of a strongly continuous semigroup 
(S(t))t >o on a Banach space X is the operator 

Ax := £r(0) = lim —(S(h)x — x) 
h\. o h 

defined for every x in its domain 

V(A) := {x G X : t i-> £ x (i) is right differentiable in t at t = 0}. 

Note that if £ x (t) is right differentiable in t at t = 0, it is also differentiable at t, for any 
t > 0. 

Theorem A.l (Hille-Yosida theorem). Let (A,T>(A)) be a positive operator on a Hilbert 
space H such that tv + A is surjective for some tv > 0. Then —A generates a contraction 
semigroup. 

We recall that on a Hilbert space, the linear operator A is called positive if (Ax, x) > 0 
for all x € T>(A). 

Theorem A.2 (Uniqueness Theorem). Let (A,T>(A)) be a closed, densely defined operator 
on a Banach space X, and let Y be a subspace of X which is continuously embedded in X 
(in symbols: Y ^ X). The part of A in Y is the operator A\ defined by A\y := Ay with 
domain T)(A\) = {y £ T>(A) n Y : Ay £ Y}. Suppose that (A,V(A)) generates a strongly 
continuous semigroup (S(t))t >o on X, and (A\,V(A\)) also generates a strongly continuous 
semigroup (R(t))t >o on Y. Then 

(A.2) S(t)y = R(t)y 

holds for all y £ Y and t > 0. Furthermore, the subspace Y is A-admissible, i.e. Y is 
an invariant subspace of S(t), t > 0, and the restriction of S(t) to Y, which is R(t), is a 
strongly continuous semigroup on Y. 

Proof. The identity (A.2) immediately follows from the uniqueness of the following Cauchy 
problem 

f ii(t) = Au(t), V t > 0, 

\«(o) = y e Y-, 

and that Y is A-admissible follows from (A.2) and the assumption. 

□ 
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Theorem A.3 (Bounded Perturbation Theorem I). Let (A,T>(A)) be the infinitesimal gen¬ 
erator of a strongly continuous semigroup (S(t))t >o on a Banach space X satisfying 

||5(i)|| < M 0 e ut , Vt > 0, 

where u € R, Mq > 1. If B € C(X), then C := A + B, with T>(C) := P(A) generates a 
strongly continuous semigroup (R(t )) t >o satisfying 

||R(f)|| < AL 0 e (w+Mo ll B ll )t , Vt > 0. 

The following results are taken from [Paz83, Chapter 5]. 

Definition A.4. A two parameter family of bounded linear operators Wft, s), 0 < s < f < 
T, on a Banach space X is called an evolution system if the two following conditions are 
satisfied: 

i) W(s 1 s) = /, Wft, r)W(r, s ) = Wft, s) for 0 <s<r<t<T; 

ii) ft, s) i->- W (t, s) is strongly continuous for 0 < s < t < T. 

Let I be the interval [0, T], and (Aft), T>(A(t))t£i be the family of infinitesimal generators 
of strongly continuous semigroups St(s), s > 0, on a Banach space X, and let Y be a Banach 
space which is densely and continuously embedded into X. The following stability definition 
appeared in [Kat70, Kat73]: 

Definition A.5 (Kato-stability and Kato-condition). We say that the family {A(t)}tg/ is 
Kato-stable, if there exist constants M > 1 and such that 

|| n n =ie s J A(t,)|| < Me^"=i s >, 

holds for every time-ordered sequence (ti,--- ,t n ) in I and Sj > 0; and we say that the 
family {A(f)} ie / satisfies the Kato-condition on I if the following conditions are satisfied: 

i) {A(t)}t£i is Kato-stable in X. 

ii) Y is A(t)-admissible for all t £ /, and the family {A(f)|y} of the part of A(t) in Y is 
Kato-stable in Y. 

in) Y C V(A(t)) holds for all t € I, and for all t £ /, Aft) is a bounded operator from Y 
into X and the mapping t > Aft) is continuous in the C(Y,X) norm || • || y x . 

Remark A.l. If for all t € I, Aft) is the infinitesimal generator of a quasi-contraction 
(see Definition A. 2) semigroup Rt(s) satisfying ||i2t(s)|| < e us , then the family {A(f)} te / is 
clearly Kato-stable with constants M = 1 and u. 

Theorem A.4 (Bounded Perturbation Theorem II). Let {A(t)}t£i be a Kato-stable family 
of infinitesimal generators with constants M and u. Let {B(t)}t£i be bounded linear oper¬ 
ators on X. If ||_£>(t)|| < K for all t 6 I, then {Aft) + B(t)} t€ i is a Kato-stable family of 
infinitesimal generators with constants M and u + KM. 

Theorem A.5. If the family {A(t)}t g / satisfies the Kato-condition (see Definition A.5) 
then there exists a unique evolution system Wft, s), 0 < s < t < T, in X satisfying 

(Ei) ||W(t,s)|| < V0 < s < t < T, 

(E 2 ) d+Wft,s)v\ t=s 


A(s)v, \/v £Y,0 < s <T, 
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(£ 3 ) d 8 W(t,s)v = —W(t, s)A(s)v, VveY, 0 <s<t<T 

Here the derivaties in (E 2 ) and (E 3 ) are in the strong sense in X and in (E 2 ) the derivative 
is from the right. 

Theorem A. 6 . Let the family {A(t)}t£i satisfies the Kato-condition (see Definition A.5) 
and let W(t, s),0 < s < t < T be the evolution system given in Theorem A.5. Suppose that 
W(t,s) further satisfies both (£4) and (£5), where 

(£4) W(t, s)Y C Y, VO < a < t < T, 

(£5) (t, s ) i->- W(t, s)v is continuous in Y for all v € Y. 

Let f £ C([s, T];Y), then for every v £ Y, u(t ) = W ( t , s)v + f*W (t, r)f(r)dr is the unique 
Y-valued solution of the following initial value problem: 

I = + /(*)> V0 < s < t < T, 

[•u(s) = v; 

and u belongs to C([s, T]; Y) n C 1 ([s, T]; X). 

Remark A.2. The definition of Y-valued solution only requires that u belongs toC(\s , T];Y )n 
C 1 ((s,T]; X) (see Definition f.l in [Paz83]j, but from the proof of Theorem 5.2 in [Paz83], 
we actually have u £ C([s, T]; Y) fl C 1 ([s, T\;X). 

Appendix B. Classical lemmas 

In this appendix, we collect some essential ingredients for Sobolev spaces (see e.g. Chapter 
13 in [Tay97] or Appendix C in [BS07]). 

Lemma B.l. Assume that U is a regular open set ofW 1 , where d is the dimension of the 
space. 

i ) Consider u and v which both belong to L°°(U) fl H S (U) with s > 0. Then their product 
also belongs to H S (U ) and there exists C > 0 depending only on s and U such that 

— c dMU ( z<olMI^(w) ~i~ IML°°(«)IMI/p>(w))' 

If s > d/2, then the L°° assumption automatically follows from the Sobolev embedding, 
and we have the following estimate: 

ll W ll/! s (W) — ( ^ll' lt ll_H' s (W)ll U ll_H' s (W)‘ 

ii) Let J 7 be a C°° function on M such that £(0) = 0. Then there exists a continuous 
function C : [0, + 00 ) i-a [0, + 00 ) such that for all u £ H S (U) fl L°°(U) with s>0: 

If s > d/2, then the L°° assumption automatically follows from the Sobolev embedding, 
and if furthermore we assume that U is bounded and that u is positive away from 0, 
i.e. |ri| > eo for some positive eo, then we have 

- € H S (U), 

u 
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if we choose F to be a C°° function such that F(x) = 0 for |x| < cq/2 and F(x) = 1/x 
for |x| > e 0 . 

in) Ifk is an integer greater than d/2+1 and a is a d-tuple of length |a| € [1 ,k\, there exists 
C > 0 depending only on k and U such that for all a in H k (U) and all u £ H\ a \~ l {U), 
we have the following estimate: 

||[<9 i a ] u || L 2( W ) < ^11 a ll IIHI rH-Rw)' 
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